Mon. Not. R. Astron. Soc. 000, 1-7 (2006) Printed 5 February 2008 (MN 1*1^ style file v2.2) 

The (^^Fe+^^Ni)/^^Ni ratio as a second parameter for Type 
la supernova properties 
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ABSTRACT 



A variation of the relative content of (^^Fe+^^Ni) versus ^^Ni may be responsi- 
ble for the observed scatter of Type la Supernovae (SNe la) about a mean relation 
between their intrinsic brightness and the shape of their light curve. Synthetic light 
curves are computed of parametrised Chandrasekhar-mass explosion models of con- 
stant kinetic energy, where the ejecta are divided into an inner NSE zone, composed of 
(^''Fe-l-^^Ni) inside and ^^Ni outside, an outer zone with Intermediate Mass Elements 
and a CO zone. Both the size of the NSE zone and the fraction of (^''Fe-h'^^Ni) v. ^'^Ni 
are varied systematically. Models with the same original NSE content but different 
(^''Fe-|-^^Ni)/'^^Ni ratios reach different peak brightness but have similar light curve 
shapes. Synthetic spectra indicate that the l^-band decline rate is not affected by the 
(54pg_^58]v^--)^56]sjj j.g^^-Q^ While the ^'^Ni mass and the total NSE mass are the dominant 
parameters determining the peak luminosity and the shape of the light curve, respec- 
tively, a variation in the (^'*Fe-|-^*Ni)/^^Ni ratio, which may depend on the metallicity 
of the progenitor (Timmes, Brown & Truran 2003) is likely to account for a significant 
part of the observed scatter of local SNe la about the mean brightness-decline rate 
relation. 
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1 INTRODUCTION 

The use of Type la Supernovae (SNe la) to constrain cos- 
mological parameters relies on the calibration of their lu- 
minosity. Although SNe la are not standard candles, sev- 
eral of their observed properties can be linked to their peak 
brightness. These properties include the decline in S-band 
magnitude in the first 15 days after B maximum, Ami5(_B) 
(Phillips 1993; Phillips et al. 1999), the shape of the light 
curve (e.g. MLCS (Riess et al. 1998), "stretch" (Perlmutter 
et al. 1997)), the ratio of two Si ll lines (Nugent et al. 1995), 
the width of the nebular emission lines (Mazzali et al. 1998) , 
the Hubble type of the parent galaxy (Hamuy et al. 2000; 
Branch, Romanishin & Baron 1996; van den Bergh 1997; van 
den Bergh, Li & Filippenko 2005), the distance from the cen- 
tre of the host galaxy (Wang, Hoflich & Wheeler 1997), and 
the B — V colour ~ 12 days after B maximum (Wang et al. 
2005). 

AU these relations indicate that the properties of SNe la 
are determined by one dominant, driving parameter, which 
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should be closely related to the amount of ^^Ni synthetised 
in the explosion (Hoflich et al. 1996). However, scatter about 
the mean relation is always present. This scatter may be 
due to uncertainties in distance and reddening affecting the 
estimate of the SN brightness, but it remains even when 
distance- and reddening- free quantities are used (e.g., neb- 
ular line width vs. A7ni5(_B); Mazzali et al. 1998). Apart 
from possible uncertainties in the measurements, this scatter 
suggests that a single parameter is not sufflcient to describe 
all possible variations of the SN la phenomenon. This has 
also been conflrmed by detailed studies of individual SNe 
(Benetti et al. 2004). 



The observed scatter could be due to a number of fac- 
tors, that may in turn be related to properties of the progen- 
itor system (e.g. distribution of ^^Ni in the ejecta, mass and 
composition of the exploding white dwarf). The identifica- 
tion of a possible second parameter may help us understand 
how SNe la explode and improve the calibration of their lu- 
minosity. Here we focus on the possible role of a variation 
of the ratio R 
Timmes et al. 



=(5''Fe-h' 
(2003). 



Ni)/ Ni, following a suggestion by 
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2 LIGHT CURVE PROPERTIES AND THE 
ROLE OF THE ("^*Fe+^*Ni)/"^«Ni RATIO 

The peak brightness of a SN la is proportional to the mass 
of ^''Ni synthesized (Arnett 1982, 1996). In fact, the decay of 
^^Ni into ^®Co and ^^Fe produces 7-rays and positrons that 
deposit their energy in the SN ejecta and eventually cause 
it to shine. It also follows from simple dimensional analysis 
that the light curve width tlc depends on the ejected mass 
Mej, the kinetic energy of the explosion Ek, and the grey 



opacity to optical photons Kopt as 

1/2 , ,3/4^-1/4 



(1) 



If the white dwarf mass and the explosion energy are as- 
sumed to be constant, the shape of the light curve must 
depend primarily on the opacity. It is now understood 
(Khokhlov, Miiller & Hoflich 1993; Hoflich et al. 1996) that 
the two quantities, mass of *®Ni and opacity, are intimately 
related, so that SN la light curves behave to first order as a 
one-parameter family. 

The thermonuclear explosion that gives rise to a SN la 
(e.g. Nomoto et al. 1984) is thought to produce material in 
nuclear statistical equilibrium (NSE) in the inner parts of 
the white dwarf. Most of this should be iron-peak nuclei. In 
the innermost ~ 0.2Mq , where the density is high enough for 
weak interactions to be important, ^*Ni and ^^Fe are mainly 
produced. In mass shells between ~ 0.2Mq and ~ O.SMq, 
weak interactions operate on a timescale long compared to 
the nuclear burning timescale, so that the electron fraction 
Ye remains constant in this region during the burning phase; 
as a consequence, the final relative abundances of ^*Ni, ^^Ni, 
and ^^Fe in this layer are directly related to the value of Ye 
at the time of the explosion. The amount of ^®Ni produced 
is essential in determining the brightness of the SN. 

Further out still, the densities are too low to burn to 
NSE, and intermediate-mass elements (IME) such as Si and 
S are produced. While no 7-rays and positrons are generated 
here, burning to Si produces almost as much kinetic en- 
ergy as burning to NSE (Gamezo, Khokhlov & Gran 2005). 
Therefore this region is also important in determining the 
properties of the SN. Finally, an outer region consisting of 
unburned C and O may also exist. 

Producing a larger amount of ''^Ni leads to more heat- 
ing. Khokhlov et al. (1993) and Hoflich et al. (1996) showed 
that the opacity increases rapidly with temperature for 
T < 10* K, the effective temperature of a typical SN la near 
maximum, but is relatively insensitive to T for higher val- 
ues. Moreover, since most ^^Ni is not mixed, a larger mass 
of ^^Ni does not immediately translate into a higher tem- 
perature outside the ^^Ni zone. 

The explanation of this behaviour lies in the fact that, 
since continuum processes are depressed in the ejecta of a 
SN la because of the lack of hydrogen and helium, the opac- 
ity is dominated by line opacity of low ionisation species 
(Pauldrach et al. 1996). Combined with the large differential 
motion of the SN envelope, line absorption causes a pseudo- 
continuum opacity which is essential in SN la light curves. 
Owing to their complex atomic level structure, ions of Fe- 
group elements (e.g. Fe 11, Fe iii, Ni 11, Co 11, Ti 11, Cr 11) have 
many more active lines than IME ions (e.g. Si 11, S 11) and 
make the dominant contribution to the line opacity (Mazzali 
et al. 2001). 



Therefore, while producing more ^ Ni means producing 
more photons (i.e. a brighter maximum), the opacity and 
the shape of the light curve depend on the total amount 
of NSE material synthesised. The line opacity is in fact the 
same in different isotopes of the same element, (e.g. ^^Fe and 
^*Fe). Since brighter SNe la have broader light curves, the 
implication is that they do not only make more ^®Ni, but 
also more NSE material. This suggests that brighter SNe do 
not simply make more ''^Ni at the expense of ^*Ni or ^''Fe. 

Mazzali et al. (2001) built parametrised explosion mod- 
els where the total NSE plus IME content is constant, so that 
Ek can also be regarded as constant, but different amounts 
of ^®Ni are synthesised at the expense of IME. This was 
achieved by moving the outer boundary of the ^®Ni zone. 
Synthetic bolometric light curves of these models reproduced 
the observed MBoi(Max)-A7ni5(_B) relation (Contardo, Lei- 
bundgut & Vacca 2000). Using spectrum synthesis, Mazzali 
et al. (2001) showed that the relation between Ms (Max) 
and Ami5{B) can also be recovered. 

The role of ''^Ni and ''^Fe is only to provide an opaque 
inner zone. Since this zone is very well hidden from the outer 
layers, the effect of its presence on the early light curve is 
limited. In this paper we explore the effect of different R val- 
ues for the same total mass of NSE material. Timmes et al. 
(2003) estimate that the final ^®Ni mass fraction depends 
on the initial metallicity Z as X('^*'Ni)« 1 - 0.054^/^0, 
a relation confirmed by detailed three-dimensional explo- 
sion and nucleosynthesis calculations by Ropke et al. (2005) 
and Travaglio, Hillebrandt & Reinecke (2005) . Both Hoflich, 
Wheeler & Thielemann (1998) and Dominguez, Hoflich & 
Straniero (2001), using delayed detonation models, also sug- 
gest that R increases with the metallicity of the progenitor. 
They find that their monochromatic and bolometric light 
curves are very little affected by changes in R. We explore 
the dependence in a more systematic way, independent of 
a specific model of the explosion. Therefore, we consider a 
situation where ^*Ni and ^^Fe are present not only at the 
lowest velocities, where they dominate, but also in the *®Ni 
zone itself. 



3 LIGHT CURVE MODELLING 

Mazzali et al. (2001) built parametrised explosion models 
consisting simply of 3 zones: the innermost zone (^''Fe, ^*Ni) 
only contributes to the opacity. Outside of this, a ''^Ni zone 
produces 7-rays and makes a large contribution to the opac- 
ity. The outer zone (IME plus CO) only makes a small con- 
tribution to the opacity. They studied three representative 
explosion models. The models have the same Mcj (I.4M0) 
and Ek (~ 1-35 x 10^^ erg), but contain different amounts of 
^^Ni: 0.4, 0.6, and O.8M0, respectively, and were accordingly 
called Ni04, Ni06, and Ni08. They reproduced the observa- 
tional properties of SNe 1992A, 1994D, and 1990N, which 
span the range of spectroscopically normal SN la decline 
rates. Starting from those models, we introduced ^*Ni and 
^*Fe in the '^^Ni layer in two different amounts, 10% and 20% 
relative to ^®Ni, and with the same distribution, following 
the predictions of Timmes et al. (2003) , and computed syn- 
thetic bolometric light curves. 

These were computed with the gray Montecarlo code 
developed by Mazzali et al. (2001). The code follows the 
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propagation and deposition of the 7-ray photons and the 
positrons emitted in the decay of ^^Ni to ''^Co and hence to 
Fe. It then follows, also in Montecarlo, the diffusion of the 
optical photons until they emerge to give rise to the observ- 
able SN light. The code adopts a simple parametrisation of 
the line opacity, based on the number of active lines in dif- 
ferent elements. In particular, it separates the effect of NSE 
elements from that of IME and of the original CO of the 
white dwarf. It also takes into account the decrease of the 
opacity after maximum caused by the cooling of the ejecta, 
as described by Khokhlov et al. (1993) and Hoflich et al. 
(1996). The opacity is parametrized as follows; 



Kopt = [0.25XFe + 0.025(1 ~ Xfc)] 



[cm^g"'], (2) 



where Xpc is the mass fraction in Fe-group elements and ta 
the time since the explosion in days. The time-dependent 
term mimicks the effect of the decreasing temperature and 
is limited to a maximum value of 2. 

Since the ^^Ni mass of the modified models is smaller 
by 10% and 20%, respectively, with respect to that of the 
original ones, it can be expected that the peak brightness of 
their respective light curves will be reduced proportionally. 
However, the effect on the opacity should be very small: 
according to Eq. (2), ^^Ni and ®'*Fe contribute to the opacity 
as much as *^Ni. In practice, we tested models with different 
*^Ni masses, but with the same amount and distribution of 
opacity. The bolometric light curves thus obtained are shown 
in Figure 1, and their most important properties are listed 
in Table 1. 

For each set of light curves in Figure 1 , the original light 
curve is the brightest one. As the ®*Ni and ^''Fe content 
increases and the ^^Ni content decreases, the light curves 
become fainter as expected. However, since models with the 
same mass of NSE material have similar opacities, their light 
curves have similar shapes. In contrast, for example, the 
Ni08-20% model light curve is similar to that of the Ni06 
model in peak brightness, but the two light curves have dif- 
ferent decline rates. Although the bolometric decline rates 
are in good agreement with observed values for SNe la cover- 
ing the range of normal spectral properties (Contardo et al. 
2000), the occurrence of SNe with different R values can gen- 
erate a dispersion about the mean MBoi(Max)-A7ni5(Bol) 
relation. Note that this would not be the case if the changes 
in composition had been introduced by extending the inner 
(^*Ni + ^*Fe) zone outwards: this would lead to the aver- 
age location of ^®Ni moving outwards as its mass becomes 
smaller and hence to a smaller photon diffusion time. This 
family of light curves would therefore follow the brightness- 
decline rate relation. 



4 EFFECT ON MONOCHROMATIC DECLINE 
RATES 

Our bolometric light curves were produced using a simple 
code, whose main advantage is that it makes it easy to trace 
the effect of properties such as the ^®Ni distribution and that 
it is not subject to uncertainties arising from the use of opac- 
ity tables. Obviously, it would be very interesting to verify 
how the monochromatic light curves are affected, as these 
are directly observable. As was done in Mazzali et al. (2001), 



we checked the behaviour of the monochromatic bands com- 
puting spectra at the appropriate epochs (maximum and 15 
days later) with our Montecarlo spectrum synthesis code 
(MazzaU & Lucy 1993; Lucy 1999; Mazzah 2000). The bolo- 
metric luminosities and the photospheric velocities obtained 
from our synthetic light curves were used as input for the 
spectral calculations. 

The results for the spectra at maximum are shown in 
Table 2. Interestingly, although the luminosity of the mod- 
els with a higher ^^Ni and ^''Fe content is smaller than that 
of the original models, the colours of the spectra are essen- 
tially unchanged, i.e. both B and V change by roughly equal 
amounts. This confirms the findings of Hofiich et al. (1998) 
and Dominguez et al. (2001) and provides further support to 
the notion that the typical colour of normal SNe la at max- 
imum light is {B — V) ~ 0.0-0.1. Furthermore, the brightest 
SNe (model NiOS) are not the bluest. This is due to the 
changing ionisation conditions, with strong Fe in lines de- 
veloping in the B band. This change from Fe 11 to Fe ill 
does not affect the light curve much: opacity is shifted to 
the blue, but so is the average flux since the temperature 
is higher. Originally, the three spectra resembled those of 
SNe 1992A, 1994D, and 1990N, respectively (Mazzah et al. 
2001, Figure 5). The spectra of the modified models are very 
similar, except for an offset in flux and a slight reduction in 
temperature as the mass of ^^Ni is reduced. 

The shape of the synthetic spectra 15 days after max- 
imum depends sensitively on the choice of photospheric ve- 
locity, which is very uncertain at such advanced epochs. This 
affects in particular the B magnitudes and the B — V colour. 
However, small changes to lower velocities than those pre- 
dicted by the light curve code are sufficient to obtain spectra 
with Am,i5(B) values similar to the original models. On the 
other hand, the V magnitude is almost unaffected by the 
choice of photospheric velocity within this range. Therefore, 
the value of Ami5(V^) can be determined quite accurately. 
We computed Ami5(l^) rather than Am2o(V), the value 
measured by Hamuy et al. (1996b), because the position of 
the photosphere would be even more uncertain at a later 
epoch. The values of Ami5{V) for our model set (Table 2) 
compare well with those of observed SNe la (Hamuy et al. 
1996b, Table 1). Since Ami5{B) and Am2o(V^) are obser- 
vationally very well correlated (Hamuy et al. 1996b), the 
agreement between our computed values of Ami5(V) and 
the observed ones supports our hypothesis that Ami5{B) 
does also not change significantly with R. 

Our calculations confirm that SNe with constant total 
NSE content but (^*Ni -I- '^■*Fe)/^'^Ni ratios varying within 
the range suggested by Timmes et al. (2003) span a range of 
~ 0.25 mag in peak brightness. However, they have similar 
V- and B-band decline rates. If SNe with a variable R do ex- 
ist, they would show a scatter in the peak magnitude-decline 
rate diagram comparable to the observed one. This is illus- 
trated in Figure 2, where we plot the location of our models 
in this diagram and the spread of observed SN properties: 
our models are able to reproduce both the brightness-decline 
relation and the observed scatter. This suggests that R may 
be one of the most important parameters responsible for the 
dispersion in the brightness-decline relation. 

More accurate pre-SN calculations and detailed explo- 
sion models may be able to determine the value and range of 
the ratio (^*Ni + ^''Fe)/^®Ni depending on properties of the 
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progenitor white dwarf such as metallicty or age, (cf. Lesaf- 
fre et aL 2006). At the same time, it will be important to 
verify the observable effect the ratio has on the observable 
properties of a SN la with more accurate radiative transfer 
calculations. Quantifying the effect and correcting for this 
may ultimately lead to a reduction of the scatter in the peak 
brightness-decline rate relation. A way to separate the first 
and the second parameter for nearby SNe could be the fol- 
lowing. On the one hand, the mass of ^^Ni can be derived 
from the late light curves, when the SN brightness is more 
directly sensitive to 7-ray and positron deposition. On the 
other hand, nebular spectra could be used to obtain infor- 
mation on both the mass of ^^Ni and that of non-radioactive 
Fe-group isotopes. 
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Table 1. Parameters of the synthetic Hght curves. 



model 


M(56Ni) 
[Mq] 


tMax 

[d] 


A^{Bol)Ma>c 

[mag] 


M(Bol)i5 
[mag] 


Ami5(Bol) 
[mag] 


Ni08 


0.80 


18.5 


-19.33 


-18.40 


0.93 


Ni08-10% 


0.72 


18.5 


-19.21 


-18.28 


0.93 


Ni08-20% 


0.64 


18.5 


-19.08 


-18.10 


0.92 


Ni06 


0.60 


18.0 


-19.15 


-18.16 


1.01 


Ni06-10% 


0.54 


18.0 


-19.03 


-18.05 


0.98 


Ni06-20% 


0.48 


18.0 


-18.90 


-17.92 


0.98 


Ni04 


0.40 


18.0 


-18.87 


-17.82 


1.05 


Ni04-10% 


0.36 


18.0 


-18.76 


-17.71 


1.05 


Ni04-20% 


0.32 


18.0 


-18.63 


-17.58 


1.05 



Table 2. Properties of the spectra and inferred value of the decline rates. 



model 


M{B)Max 


{B - y)Max 


(B-V)i5 


Ami5(B) 


Ami5(y) 




[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


Ni08 


-19.22 


0.10 


0.56 


1.08 


0.62 


Ni08-10% 


-19.08 


0.10 


0.61 


1.08 


0.57 


Ni08-20% 


-18.94 


0.16 


0.65 


1.09 


0.60 


Ni06 


-19.12 


0.06 


0.70 


1.31 


0.67 


Ni06-10% 


-19.03 


0.08 


0.71 


1.30 


0.67 


Ni06-20% 


-18.85 


0.08 


0.79 


1.32 


0.61 


Ni04 


-18.88 


0.06 


0.81 


1.46 


0.71 


Ni04-10% 


-18.76 


0.06 


0.84 


1.46 


0.68 


Ni04-20% 


-18.63 


0.07 


0.81 


1.45 


0.71 
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Figure 1. Families of bolometric light curves computed for models Ni08, Ni06 and Ni04, with and without the inclusion of ®^Ni and 
^''Fc. The continuous curves are for the original models (with O.SMq, O.GMq, and O.4M0 of ^®Ni, respectively), the dotted ones are for 
compositions where 10% of ^^Ni has been replaced with ^*Ni and ^*Fe, while in the dashed ones this value is 20%. In all cases, light 
curves based on model Ni08 are the brightest, while those based on model Ni04 are the dimmest. 
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Figure 2. The position of the models on the brightness-decline rate diagram. The ridge line is the locus of the Calan-Tololo SNc (Hamuy 
et al. 1996a), shifted for Hq = 72kms~^ Mpc~^, and the dashed lines are the l-cr dispersion. 



